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ABSTRACT 


The eddy diffusion coefficient is estimated as a function of altitude, separately for the Jovian 
troposphere and mesosphere. The growth-rate and motion of particles is estimated for various 
substances: the water clouds are probably nucleated by NH,Cl, and sodium compounds are likely 
to be absent at and above the levels of the water clouds. Complex organic molecules produced by 
the Læ photolysis of methane may possibly be the absorbers in the lower mesosphere which 
account for the low reflectivity of Jupiter in the near-ultraviolet. The optical frequency chromo- 
phores are localized at or just below the Jovian tropopause. Candidate chromophore molecules 
must satisfy the condition that they are produced sufficiently rapidly that convective pyrolysis 
maintains the observed chromophore optical depth. Organic molecules and polymeric sulfur 
produced through HS photolysis at A > 2300 Å probably fail this test, even if a slow, deep 
circulation pattern, driven by latent heat, is present. The condition may be satisfied if complex 
organic chromophores are produced with high quantum yield by NH; photolysis at A < 2300 A. 
However, Jovian photoautotrophs in the upper troposphere satisfy this condition well, even with 
fast circulation, assuming only biochemical properties of comparable terrestrial organisms. Unless 
buoyancy can be achieved, a hypothetical organism drifts downward and is pyrolyzed. An 
organism in the form of a thin, gas-filled balloon can grow fast enough to replicate if (i) it can 
survive at the low mesospheric temperatures, or if (ii) photosynthesis occurs in the troposphere. 
If hypothetical organisms are capable of slow, powered locomotion and coalescence, they can 
grow large enough to achieve buoyancy. Ecological niches for sinkers, floaters, and hunters 


appear to exist in the Jovian atmosphere. 


Subject headings: planets: atmospheres — planets: Jupiter 


I. INTRODUCTION 


The contemporary Jovian atmosphere has some 
similarities to the primitive terrestrial atmosphere, 
where organic molecules were produced readily; and 
experiments in which a wide variety of organic 
molecules are produced under conditions similar to 
those of the Jovian atmosphere have been carried out 
(Sagan and Miller 1960; Sagan et al. 1967; Woeller 
and Ponnamperuma 1969; Rabinowitz et al. 1969; 
Sagan and Khare 197la, b; Khare and Sagan 1973; 
Ferris and Chen 1975). This has led to the hypothesis 
that there may be a Jovian biology (Sagan 1961; 
Shklovskii and Sagan 1966). Observations in the 
visible and in the ultraviolet have implied the presence 
of particulate matter in the Jovian atmosphere which 
may (but need not necessarily) be connected with 
questions of Jovian organic chemistry and biology. 
Both chemical and biological issues are affected by 
fluid dynamics in the Jovian atmosphere. We are here 
concerned largely with these effects. 

Modern models of the structure of Jupiter (for 
references, see, e.g., Stevenson and Salpeter 1976) 
predict that the bulk constituent, a hydrogen-helium 
mixture, is fluid throughout, with temperatures in- 
creasing inward to above 10* K. Plausible organisms 
require temperatures well below 10° K. One principal 
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problem which must be faced by such hypothetical 
organisms is descent to pyrolytic depths due either to 
the acceleration of gravity or to convective down- 
drafts. The second half of the present paper examines 
(1) whether such organisms undergo significant growth 
and replication in less than the time scale for drift to 
pyrolytic depths, or (2) whether other adaptations to 
avoid pyrolysis are hydrodynamically and biologically 
feasible. We shall see that (1) becomes easier for 
smaller organisms and (2) for larger ones. 

That plausible and internally consistent Jovian 
ecologies can be described cannot, by itself, demon- 
strate the likelihood of life on Jupiter. A tenable 
argument for life on Jupiter must also show that the 
origin of life on that planet was possible. Since many 
steps, particularly the later ones, which led to the 
origin of life on Earth are only incompletely under- 
stood (see, e.g., Miller and Orgel 1973), a thorough 
discussion of the origin of life on Jupiter cannot be 
undertaken at the present time. Most skepticism on 
this subject seems to arise from the contention that 
synthesized molecules will be carried convectively to 
pyrolytic depths before the origin of life can occur. 
The unspoken premise in such an argument is that 
pyrolysis would not have occurred on the primitive 
Earth. However, the measured Arrhenius rate con- 
stants imply, through a regression analysis, that many 
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essential organic molecules, for example a number of 
biological amino acids, suffer thermal degradation at 
typical terrestrial surface temperatures in geologically 
short time periods. The very limited available experi- 
mental data show, for example, that the half-life at 
30° C of the amino acid serine is ~10* yrs and of 
threonine is < 10° yrs; while at the same temperature 
the half-life of the simplest amino acid, alanine, is 
>10?° yrs (Vallentyne 1964). Because of the ex- 
ponential Arrhenius kinetics, the highest temperatures 
on a temperature-variable planet dominate the net 
molecular degradations; also there is evidence that 
aqueous solution and alkaline pH’s—both expected on 
the primitive Earth—enhance the decomposition rates 
(Vallentyne 1964). It does not seem at all improbable 
that some essential precursor molecules for the origin 
of life had half-lives on the primitive Earth of 10? to 
10° yrs, not greatly dissimilar to the lifetimes against 
convective pyrolysis of organic molecules on con- 
temporary Jupiter, as calculated below. 

Many workers on the origin of life have suggested 
that the rate-limiting steps were stochastic in nature. 
The volume available for natural molecular experi- 
ments on Jupiter is approximately 10° times larger 
than on Earth. If phase interfaces play a critical role 
in the origin of life (see, e.g., Bernal 1973) Jupiter, 
with a great concentration of solid and liquid at- 
mospheric aerosol particles, would be particularly 
favorable. 

The origin of life on Earth took a period <0.1 the 
lifetime of the Earth, and possibly much less (Sagan 
1975). On Jupiter, we have much longer periods of 
time available for chemical experiments, a larger 
interaction volume, and a much larger area for phase 
interface chemistry. On the other hand, the thermal 
degradation time scales are shorter for many mole- 
cules, and there is a much greater excess of molecular 
hydrogen than on the primitive Earth, where H, was 
rapidly depleted during planetary formation and in 
subsequent Jeans escape. The effect of molecular 
hydrogen on prebiological organic chemistry is in- 
adequately understood, but is briefly discussed below. 
Its effect on quenching gas-phase organic reactions is 
less than linear, at least for some reactions. We con- 
clude that the absolute a priori likelihood of life on 
Jupiter, as on Earth, cannot be estimated reliably at 
the present time, but, as on Earth, can hardly be ruled 
out. The relative a priori probabilities of the origin of 
life, in a comparison of Jupiter and Earth, depend on 
which reactions and molecular precursors are rate- 
limiting. If, for example, the alanine abundance is 
rate-limiting, the Earth is strongly favored; if serine is 
rate-limiting, Jupiter may be favored. It is also possible 
that there exists a variety of different paths to the origin 
of life, for some of which the Earth is favored and for 
others of which Jupiter is favored. 

We devote the remainder of the present paper to a 
characterization of the relevant physical environment 
of Jupiter; to a discussion of the chromophores 
responsible for the observed coloration of the planet; 
and to an analysis of some permissible ecological 
niches of hypothetical Jovian organisms. 
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Values for the eddy diffusion coefficient K play an 
important role in such arguments. Estimates are given 
in § II, separately for the convective troposphere and 
for the more stable mesosphere. Analytic approxima- 
tions are also given there for abundance as a function 
of height for particles which are produced at some 
upper level and are transported downward toward 
pyrolysis. 

We are concerned with two sorts of particulate 
matter: particles that condense out under conditions 
of thermal equilibrium, and particles that require an 
investment of free energy for their production. 
Equilibrium condensation has been studied previously 
for H2O, NH,SH, and NH3; we investigate it further 
in § III for a number of substances with higher con- 
densation temperatures. We estimate growth rates for 
these condensates and discuss whether NaCl or other 
particles can be present at the level of the water clouds 
to act as nucleation agents. The production of 
chromophores and of nonequilibrium organic mole- 
cules is treated in § IV. We obtain a lower limit to the 
production rate of ethane from methane photolysis by 
solar La, and treat organic chromophores produced by 
various wavelengths of light, as well as polymeric 
sulfur chromophores. 

The motion of hypothetical organisms is discussed 
in § V, and the diffusion of metabolites and the 
consequent growth of organisms in § VI. We consider 
four different cases, depending on whether organisms 
can or cannot survive at the low temperatures of the 
mesosphere and whether they can or cannot utilize 
visible light for photosynthesis. In § VII we consider 
the more rapid growth attainable if organisms can 
control their motions. 


II. EDDY DIFFUSION IN THE JOVIAN ATMOSPHERE 


The total infrared flux emerging from Jupiter 
(Aumann, Gillespie, and Low 1969; Ingersoll et al. 
1975) is Faos = oTa © 1.5 x 10* ergs cm~? s7? 
(Tere © 127 K), with about half of this flux coming 
from the deep interior (Fins % 0.8 x 10* ergs cm~? 
s~1). The radiative opacity increases fairly rapidly 
with increasing pressure P, and radiative-convective 
model calculations predict a fairly rapid change across 
the boundary level defined here as the tropopause. 
Above the tropopause we have radiative (subadiabatic) 
heat transport in the mesosphere; below we have con- 
vective heat transport (with the temperature gradient 
only slightly superadiabatic) in the troposphere. 
Model atmospheres (Trafton and Munch 1969; 
Hogan, Rasool, and Encrenaz 1969; Wallace, Prather, 
and Belton 1974) and spectroscopic observations 
(Gillett, Low, and Stein 1969; Aitken and Jones 1972; 
Lacy et al. 1975) suggest a temperature T x 175 K at 
the pressure level of P = 1 bar and place the tropo- 
pause somewhere slightly above this. We shall adopt 
the level (T ~ 145, P x 0.5 bar) for the tropopause. 
These numbers are consistent with Pioneer 10 infrared 
radiometer data (Orton 1975a) and with ground-based 
infrared radiometric measures (Orton 1975b) and 
possibly also with Pioneer 10 and 11 S-band occulta- 
tion data (Kliore and Woiceshyn 1976). We assume a 
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helium abundance of about 10% by number, consistent 
with separate measurements from the occultation of 
Beta Scorpii (Elliot et al. 1974); from Pioneer 10 and 11 
ultraviolet photometry (Carlson and Judge 1974); and 
from the pressure broadening of the J = 0 and J = 1 
hydrogen lines (Houck et al. 1975). Assuming a helium 
abundance of about 10% (by number), it then follows 
that y/(y — 1) © 3, where y is the ratio of specific 
heats, and that the mean molecular weight p ~ 2.2. 
For pressure P, pressure scale height H = P(dP/dz)-1, 


density p, and kinematic viscosity v in the troposphere ' 


we adopt 


175 K 


a EEL 70.5 cm? s~! (1) 
V ~ 175K . . 


In the mesosphere, the temperature gradient decreases 
rapidly (for P < 0.2 bar, say), and the minimum tem- 
perature is about 120 K followed by a temperature 
inversion whose structure is not yet well known 
(possibly Tmax % 150 K at P ~ 10-~°? bar). 

For a general convective circulation pattern the 
vertical convective heat flux F, is given by 


F, = c,pATw, , (2) 


where c, is the specific heat at constant pressure, AT is 
the typical vertical excess in the temperature drop, and 
w, the typical vertical convection speed. For vertical 
distances larger than the vertical extent / of the cir- 
culation pattern, one has an eddy diffusion coefficient 
K, = w-l; and we shall denote by t, the eddy diffusion 
time, H?/K,, for diffusion through a single pressure 
scale height H. In the simplest form of Prandtl mixing- 
length theory (see, e.g., Schwarzschild 1958), the 
circulation is characterized by two parameters, the 
mixing length / and the superadiabaticity with w,? and 
AT/T both proportional to the superadiabaticity. We 
replace F, by 1.4 x 10* ergs cm~? s~?, a typical value 
for the Jovian troposphere, employ equation (1), and 
write K., We, and t, in the forms 


K, = &(7/175 K)" x 107 cm?s7?, 
We = &/4(T/175 K)? x 76 cms7?, 
te = €-*(T/175 K)? x 7.5 x 10°s. (3) 
In the notation of mixing-length theory our parameter 
é equals (//0.05H)*’*. The Jovian tropospheric lapse 
rate is then very close to adiabatic (A7/T ~ 107°), as 


is true by the same argument (Sagan 1960) for all 
other planetary tropospheres. 
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In the simplest version of mixing-length theory 
l x Hand £é ~ 100. The actual circulation patterns in 
the troposphere are likely to be much more complex 
(Hunten 1975), In particular, one should consider the 
effects of rotation at great depths and of latent heat 
release in the cloud layers near the tropopause 
(Gierasch 1973). For a constant value of the heat flux 
F, in equation (2), these effects usually increase AT and 
decrease w, considerably, leading to smaller values of 
K, and £ in equation (3). These considerations suggest, 
as a simple working hypothesis, € = 1, which should 
be accurate to within a factor of about 200. 

Gierasch (1976) has suggested a specific large-scale 
circulation pattern which extends all the way from the 
bottom of the water-cloud layer (T ~ 280 K) to the 
tropopause. This pattern is driven by the latent heat 
of water which leads to a relatively large value of 
AT/T ~ 10-?. This corresponds to é ~ 0.01 with 
time scales for complete turnover of the region of up 
to a few years, even though some cloud patterns 
(above the water clouds) are observed to change in 
shorter times (Peek 1958). On Gierasch’s model, the 
eddy diffusion coefficient increases by a very large 
factor just below the cloud bottom to € ~ 10. 

In the mesosphere the bulk of the heat flux F,,, is 
carried by radiation, but some fraction of the flux 
could drive a circulation pattern working against a 
subadiabatic temperature gradient. The vertical speed 
Wm Of such a circulation pattern must satisfy an in- 
equality, based on the fact that the adiabatic cooling 
rate cannot exceed the total available heating rate: 


Hpw,,[c,(a4T/dz) + g] < Faot- (4) 


Except in the immediate vicinity of the tropopause, the 
temperature gradient in equation (4) can be neglected 
for the troposphere. Putting T ~ 150 K and assuming 
that the vertical length scale for the circulation pattern 
is of order H, the eddy diffusion coefficient Km is of 
order w,,H, so that 


Wm < (1 bar/P)0.016 cm s~?; 
Km < (1 bar/P)7 x 10*cm?s7?. (5) 


Thus typical velocities just below the tropopause are 
many thousands of times larger than just above the 
tropopause. Analysis of the observed Le albedo of the 
upper mesosphere (Wallace and Hunten 1973; Carlson 
and Judge 1974) indicates K,, ~ 3 x 108+! cm?s~? at 
levels where P ~ 3 x 10-7 bar; analysis of B Sco 
occultation measurements (Sagan et al. 1974) indicates 
Km > 108 cm? s~? at similar levels. It is theoretically 
plausible (French and Gierasch 1974) that gravity-wave 
propagation upward in the mesosphere contributes to 
the eddy diffusion. For a single wave mode, K,, has a 
step-function behavior, but overall it seems plausible 
to assume some constant fraction of the maximum 
expression in equation (5). The value given by the Le 
albedo then suggests 


Km ~ (1 bar/P)100 cm? s~}, (6) 
uncertain by factors of order 10 or 100. Note that 
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Fic. 1.—A schematic plot of temperature T and eddy diffusion coefficient K versus pressure P for the Jovian atmosphere. 
The dashed portion of the curve for K corresponds to Gierasch’s (1976) model for a deep circulation pattern driven by latent 
heat effects. Methane photolysis extends upward from the indicated level. 


Km œ e?/#, Temperature T and K(z) are plotted against 
pressure, schematically, in Figure 1. (The dashed 
curve for K corresponds to Gierasch’s 1976 model in 
which the latent heat of water condensation is taken 
into account.) 

We shall need continuity equations for nz), the 
number density as a function of height z of some 
complex molecule (of molecular weight u, > 2.2) or of 
some particle. Let w,(z) be the steady-state vertical 
drift velocity of the molecule or droplet falling under 
gravity through an atmosphere with pressure scale 
height H, eddy diffusion coefficient K, and mean 
molecular weight 2.2. In standard notation, the con- 
tinuity equation (see, e.g., Strobel 1973) reads 


— 


j’ 


d K D; 
(K + D) = + (5+ y= 


where D, is the effective molecular (or particle) 
diffusion coefficient and H, = A(2.2/u;) « H is the 
effective scale height for the jth particles alone. 9, is 
the vertical flux of particles of species j. H, and the 
drift velocity w4 for free fall under gravity are related 
by D, = waH;. We shall be interested only in levels 
well below the turbopause where D, « K and give 
only a simplified version of the standard treatment 
(Strobel 1973; Prinn 1973). The steady-state continuity 
equation then reduces to 


Kt 5 AP + woho- - 0G). o 


In § IV we shall apply equation (7) to the number 
density n(z) of organic molecules or other chromophore 
particles for which the drift velocity can be neglected 
(w; « K/H). Assume further that the particles are 
mainly produced at some higher level z % Zo > 0 and 
destroyed at some lower level z = 0, so that the flux ® 
equals the total production rate (cm~?s~*) for 
0<z<Z, and ® = 0 outside that range. If the 
production level Zo is in the troposphere, the eddy 
diffusion coefficient is given by K, in equation (3), 
which varies quite slowly over a scale height H. An 
approximate solution of equation (7), which neglects 
terms in dK,/dz, is then 


n(z)  (PA/K,)[1 — exp(—2/H)] for0 <z< Zo, 
x (DHIK) exp [—(z — Z,)/H] for H «K Zo <z. 
(8) 


Thus particle number densities decline exponentially 
with altitude above the production level and decrease 
more slowly with depth below the production level. 
Mathematically, the solution in equation (8) depends 
on the boundary condition that n(z) = 0 for z < 0 
and that © is already constant for infinitesimally small 
positive z. Physically, the requirement is merely that 
the destruction rate change from very small values (at 
positive z) to very large values (at negative z) over less 
than a scale height. This should be a good approxima- 
tion since pyrolysis rates are very temperature- 
sensitive. (But note that the level z = 0 refers to 
different altitudes for different molecular species.) 
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The situation is more complicated if the production 
level Zo is in the mesosphere and the destruction level 
(z = 0) in the troposphere, with z = Z, < Zo repre- 
senting the tropopause. However, neglecting again 
derivatives of K = K, (for z < Z,) and assuming that 
K = Kn (for Z, < z) is proportional to e”¥ as in 
equation (6), the continuity equation (7) can again be 
solved analytically. With the fact that Z, > H and 
KZ») > Kn(Z,) for cases of practical interest, the 
solution can be approximated by 


n(z) x (®H/K.)[1 — exp(—z/H)] for0 <z<Zp, 


-z+ HEA 


ae 
~ KZ) 
for Zp < Z < Zo, 
x [®(Zo — Z,)/Kn(Zo)] exp [- (z — Zo/H)] 
for Zo < z. (9) 


. Thus particle number densities increase slowly with 
altitude above the destruction level until the tropo- 
pause, decline slower than exponentially (ocze~?/*) 
between the tropopause and the production level, and 
decline exponentially above the production level. For 
mesospheric production and tropospheric destruction, 
we therefore expect maximum number densities just 
above the tropopause. Again, no particles exist at 
z < 0. Note that the mixing ratio (ocn/p,,,) remains 
constant above the production level whether this level 
is in the mesosphere or troposphere. In reality, this 
mixing ratio would decrease slightly with increasing 
height z if some additional destruction takes place 
above the production level or if the drift velocity wg is 
not entirely negligible. 


HI. THERMAL CONDENSATION IN DEEP CLOUD LAYERS 


The thermal condensation of droplets to form cloud 
layers in the Jovian atmosphere has been discussed 
(Weidenschilling and Lewis 1973; Stauffer and Kiang 


' 1974) only for the upper levels of the troposphere 


(T < 300 K), where direct observations are possible. 
Molecules involving H, C, N, O, P, and S have been 
considered, which lead especially to condensation of 
H:O (T < 275 K), NH,SH (T < 220K), and NH; 
(T <= 160 K). The condensation of droplets of more 
refractory materials in deeper layers (higher tempera- 
ture, T) is also of some interest, especially for com- 
pounds containing Na and/or Cl. If such small grains 
or droplets can be convected to the upper troposphere, 
they can act as seed nuclei for the condensation of 
H,O drops (at T ~ 280 K). Their abundance is also 
relevant for any hypothetical Jovian biology and for 
the question of the replication of terrestrial microbial 
contaminants in the Jovian clouds. The following 
crude estimates and conjectures should serve at least 
for purposes of orientation. 

In the very deep interior of Jupiter where the density 
exceeds 0.2 gm cm~? (P > 10° bar, T > 10* K), in- 
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complete solubility of various substances in hydrogen 
could lead to substantial chemical separation. This 
possibility, discussed so far only for He (Salpeter 
1973; Smoluchowski 1973), should be investigated 
further; however, we assume provisionally cosmic 
abundances of all elements (Allen 1973) at P % 10°-5 
bar (T ~ 2500 K), say. At this level practically all 
substances are gaseous. Refractory molecules begin to 
condense at slightly lower temperatures. 

For each material, with f the cosmic abundance by 
number (mixing ratio) of the least abundant element in 
the compound, we must evaluate a number of 
quantities: 

1. The equilibrium condensation temperature Teon 
appropriate to the Jovian troposphere, i.e., that tem- 
perature at which the saturation vapor pressure equals 
f times the pressure in equation (1). Values for Toon 
have already been given by Lewis (1969b). 

2. Assuming the absence of any external seed 
nuclei, we calculate the amount of supercooling ATer 
required from homogeneous nucleation theory. For 
this purpose we need a parameter 7, defined in equation 
(7) of Salpeter (1974), which is essentially the number 
of times appropriate molecules of a particular species 
stick to a given surface site on a droplet during a 
typical flow time. In that definition for 7 we use 
s ~ 0.5 for the surface sticking coefficient and equa- 
tion (1) for the pressure, and replace the flow velocity 
v by the typical convection speed w, in equation (3), 
which gives n © f(Toon/175 K)*? x 3 x 10**. Values 
of f and T,,, (taken from JANAF 1971) are given in 
Table 1 for a number of compounds, and 7 exceeds 
4 x 108 in all cases. The ratio AT,,/Toon is given by 
equation (9) of Salpeter (1974), but multiplied by a 
factor 6712 (correcting an error kindly pointed out by 
B. Draine). This ratio is fairly small in the cases con- 
sidered and not of direct importance for our crude 
estimates. But we do need the final radius a, of typical 
drops which would result at the end of purely homo- 
geneous nucleation. 

At first sight one would expect a, to be of order 
times the intermolecular spacing in the liquid, but 
there are a number of correction factors. The most 
important correction is due to the fact that the drop 
radius a becomes large compared with the mean free 
path for a molecule in the gas, Im ~ (T/175 K)~? x 
6 x 10-5 cm. The droplet’s growth is controlled by 
the diffusion speed rather than the gas-kinetic speed, 
and the rate is multiplied by a correction factor 
(1,,/a) « 1. Including some further logarithmic factors, 
we find, as a very rough approximation, 


as ~ (175 K/T)y/7(logig n) x 10-& cm. (10) 


Our values for a; are given in Table 1 and have to be 
compared with a number of other length scales. 

3. Let ag, be a critical radius for a spherical droplet 
such that the hydrodynamic drag coefficient Cp is 
unity. For a « dp, the Reynolds number is small and 
the drag coefficient is given by the Stokes formula 
Cp = (12v/aw,), where wg is the Stokes expression for 
the terminal vertical drift velocity of the sphere. Using 
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TABLE 1 


EQUILIBRIUM CONDENSATION TEMPERATURE, Toon, AND ESTIMATES OF THE RELEVANT MIXING RATIO, f, THE RADIUS a;, THE 
PARAMETER ‘7jyvap, AND THE MINIMUM MIXING RATIO fm FOR A NUMBER OF SUBSTANCES ABOVE THEIR CONDENSATION LEVELS 


Parameter Fe Meg2SiO, 
logro f. wee eee cece eee 4.1 4.6 
con/ LOO K....... 2. eee eee 25 23 
ayf10-F com. wwe eee 120 60 
Trap Ce ey 60 20 
— logio fine sce cece cenececs 7.9 7.7 
equation (1), we have 
äre = 0.007 om (2 ZX) oe >, 
Par T 
2 a°g par 
wala) = -—- = 
a( ) 9 v 7 
_ T \ —1/2 Pa a 2 
~ If = Far); ; 
ee (T5 z) ( 2 I(T =a) 
(11) 


where par is the internal density (in gm cm~*) of a 
liquid drop. Another critical length is a,, the droplet 
radius for which the drift velocity w4 equals the typical 
convection speed w, in equation (3). For the cases of 
interest (T < 2500 K) we have a4 < are, so that 
equation (11) holds. With pa %2 and éz 1 in 
equation (3), we have 


aa X 0.002 cm(7/175 K)~122., (12) 


Large drops with a > a, “‘rain out” rapidly, whereas 
small droplets with a « a, are carried along by the 
convection pattern and can move up to higher levels. 

For all cases in Table 1 (except possibly for NHC) 
we have a; > az, so that homogeneous nucleation 


leads to large enough drops for most of the matter to 


be rained out. However, if a very large number of 
nucleation seeds are present originally, a large number 
of much smaller droplets condense out, and the drop- 
lets must grow before they can fall out as rain. Let 
Geo, be that radius for which Cp = pa;/p,: 


— 2/3 
dooa ~ (4) (asx) x 3 x 10-*cm. (13) 


Once @ È aa, droplets usually can grow very quickly 
by coalescence, but for a < doo, growth must occur by 
the slower process of smaller drops evaporating and 
the excess vapor diffusing to larger drops (Langmuir 
1948; Mason 1972; Rossow 1976). An important 
parameter nyap is the growth rate from this process 
(from radius 0.5đ.oa tO Ago, SAY) expressed in units of 
H?/K,, a typical convection time. With an uncertainty 
of at least a factor of 10, we have 


T 7/16 ‘ 2 
vap ~ 4 X 1047( 75) (%) i (14) 


SiO2 NazSi205 NaCl NH.Cl (NaS) 
4.4 5.7 6.3 6.3 (5.7) 
16 10 4 10 
20 4 2 8 
1 0.1 0.03 0.4 
7.2 7.8 7.0 6.4 7.0 


If vap > 1 for a particular material, most of it will 
form into large “raindrops” in a cloud layer at levels 
close to the condensation temperature Ton, and only 
small amounts of the material will be found at higher, 
cooler levels. If nyap « 1 one might expect only small 
droplets of the material to form, so that wz « w, and 
upward convection can maintain an appreciable mix- 
ing ratio of the droplets at higher levels. In practice, 
however, various complications usually tend toward 
more rapid formation of larger droplets (Rossow 
1976). For nyap ~ 1071 or 1072, the circumstances 
are probably mixed: some rain near the 7,,, level, 
some small droplets a few scale heights higher up, but 
with a mixing ratio which decreases slightly with in- 
creasing height. Values for nyap from equation (14) are 
given in Table 1. Except possibly for NH,Cl, nyap > - 
10-2 as well as a, > ag, and most of the material rains 
out and cannot move up to higher levels. However, in 
upward-moving convection cells the saturation vapor 
pressure continues to fall; droplets then continue to 
condense out, but with decreasing values of the mixing 
ratio f and hence of the parameter 7. Let fm be that 
(decreased) value of the mixing ratio, for which the 
droplet size a; (due to homogeneous nucleation) in 
equation (10) equals az. Values of fa are also given in 
Table 1. The mixing ratio will not drop appreciably 
below fm in upward moving convection cells, since 
droplets are now too small to rain out rapidly enough 
(and further growth of droplets is also too slow). 

For the cases with 7T,., > 1600 K the radii satisfy 
a; > Aeon, COalescence is rapid, and large drops are 
rained out immediately. At the high pressures in 
Jupiter’s interior, liquid Mg.SiO, (as well as Ca and 
Al silicates and metallic Fe) condenses before MgSiO, 
and before SiO,. When liquid SiO, begins condensing 
near the 1800 K level, most of the metallic elements 
have already been rained out, but sodium has not. For 
pure Na,SiO; our estimates indicate a slightly lower 
condensation temperature than for SiO, which 
introduces an uncertainty: It is conceivable that all the 
silicon has been rained out before sodium condenses, 
in which case gaseous NaCl and NaOH (since the 
elemental abundance of Na exceeds that of Cl) survive 
to higher, cooler levels. It is more likely that much of 
the sodium condenses out between 1800 K and 1600 K 
in the form of some sodium silicate-silica solution, in 
which case gaseous HCl and some NaCl survive to 
higher levels. In either case, more of the surviving 
NaCl condenses out near 1000 K; in the more likely 
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case, NH,Cl condenses near 400 K; in the less likely 
case, NaS condenses near 1000 K. 

We have some tentative implications for levels with 
T < 700 K, where Jovian biology (or at least organic 
chemistry) might be possible: Sodium (or other alkali 
metals and alkaline earths) and chlorine (or other 
halogens) cannot be present simultaneously. However, 
either sodium sulfide or (more likely) ammonium 
chloride is present. For the more likely case, some 
fraction of the NH,Cl crystals are convected upward 
and act as seed nuclei for the condensation of water 
droplets, so that essentially no supercooling is present 
in the base level of the water clouds (Stauffer and 
Kiang 1974). For the less likely case of Na.S crystals, 
convection over a larger height difference is required 
for nucleating the water, and f,, is slightly smaller. 
However, the chromophores to be discussed in § IV, 
as well as NaS and NH,Cl, could probably also act as 
seed nuclei, and it seems generally safe to assume no 
supersaturation for water. We estimate below the 
amount of exogenous sodium which may arrive at the 
cloud level from micrometeoritic infall. 

Since both sulfur and sodium have been detected in 
toroidal nebulae about Jupiter—probably arising 
from the surface of Io—these clouds are a potential 
source of S and Na in the clouds of Jupiter. From 
observations of forbidden S 11 circumjovian line 
emission, Brown (1976) calculates a total sulfur flux 
from the surface of Io of 10®-6*°-7 cm~? s+, equivalent 
to roughly 10* cm~?s~? at the Jovian clouds. This is 
some 8 orders of magnitude less than that expected 
from H.S photodissociation (see § IV), and therefore 
is entirely negligible. The cosmic abundance of Na is 
at least an order of magnitude less than that of S, and 
so only very small quantities of sodium could be 
present at the Jovian clouds due to leakage from the 
Io Na toroid. 

Approximate influxes of micrometeoritic sodium 
and sulfur, assuming the fluxes are the same on Jupiter 
as on Earth, and that micrometeorites are primarily 
cometary ices with cosmic abundances of other atoms, 
are ~10° cm~? s~! for Na, and ~10® cm~? s-t for S. 
Thus the micrometeoritic infall rate of these atoms is 
likely to exceed that from the Io-associated toroidal 
clouds. In the upper troposphere with T ~ 175K, 
equations (3) and (9) then give n ~ 10* Na atoms 
cm~%, corresponding to a mixing ratio f < 10715, and 
even smaller at the water cloud levels. These values are 
small compared with the values of fa in Table 1. 


IV. PHOTOPRODUCTION OF MOLECULES AND 
CHROMOPHORES 


a) Far-Ultraviolet 


At thermal equilibrium methane is by far the most 
abundant carbon compound in the Jovian atmosphere, 
where it remains in-gaseous form. Methane absorbs 
UV photons at all wavelengths shorter than 1600 A, 
and the total flux of such photons at the top of the 
Jovian atmosphere (mainly from solar La at A= 
1216 A) is ~5 x 10° photons cm~?s~! (global di- 
urnal average). Very little of the Le photons are 
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scattered by neutral hydrogen (Carlson and Judge 
1974), and about one-third of the photons process 
CH, into other hydrocarbons (Strobel 1973a, 1975). 
The dominant product is ethane with a predicted 
production rate of R ~ 1 x 10°C,H,cm~?s7! ~ 
5 x 10714 gm cm~% s-!, The ethane produced is 
rather stable and is destroyed mainly by eddy diffusion 
into the troposphere, followed by pyrolysis in deeper, 
hotter layers. For an assumed mesosphere eddy 
diffusion coefficient K,,(z), accurate numerical calcula- 
tions for the number density n(z) of CHs can be 
carried out (Strobel 1975). However, equation (9) is 
sufficient for purposes of orientation (the largest un- 
certainty comes from uncertainties in K,,). Assuming 
the tropopause level Z, to occur at P ~ 0.5 bar (with 
P ~ 107° to 10-77 bar at the production level Z,), and 
with K,,(z) given by equation (6), the maximum ethane 
number density occurs one scale height H above the 
tropopause: 


Nmax © DA/K(Z, + H) 
=~ 5 x 10 cm-? 25 x 10-"g,, (15) 


where K(Z, + H) is the eddy diffusion coefficient 
evaluated at altitude Z, + H. 

The maximum value of the mixing ratio occurs just 
below the production level, and the total predicted 
column density N of ethane molecules is 


N X ®H?/Kpy(Zp) 
z 5 x 10° cm~? x 3 x 1078? gcem-?. (16) 


Ethane (and smaller amounts of acetylene) have 
indeed been detected in the Jovian mesosphere by 
means of the emission spectrum at Az 10yum 
(Ridgway 1974). High-resolution spectra are not yet 
available, and the column density N of C,H, from the 
present spectra depends strongly on the distribution of 
gas temperatures, which is not well determined by 
low-resolution spectra: If the temperature in the 
mesospheric inversion layer turns out to have a 
maximum value of 150 K, the observations indicate 
(Ridgway 1974) ~0.02 g cm~? of ethane in this high- 
temperature region. This would be in very serious 
conflict with the predictions, since only a small fraction 
of the theoretical column density in equation (16) 
refers to the hotter region many scale heights above the 
tropopause. Presumably one is dealing with some 
combination of (a) higher mesospheric temperatures, 
(b) smaller diffusion rates than in equation (6), or (c) 
some additional source of ethane raising the produc- 
tion rate R above 10° cm~? s71, 

Production rates from the same Le radiation in the 
mesosphere for larger, more complex organic mole- 
cules have been estimated (Strobel 1975) to be factors 
of ~10 or 100 (by number) less than ethane, an 
estimate consistent with the data presented by Noyes 
and Leighton (1941), who remark that the ultraviolet 
irradiation of ethane leads to “‘a surprising variety of 
products,’ including unsaturated polymeric hydro- 
carbons. Assuming one-tenth the rate (by mass), the 
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mass-production rate Rog is ~5 x 10715 g cm~? s71, 
and equation (16) would give a column density for 
complex organics of Norg ~ 3 x 1074 g cm~?. 
Complex organic molecules may or may not be the 
mesospheric chromophores which have been inferred 
from the UV albedos: The measured albedo of 
Jupiter (Wallace et al. 1972; Savage and Caldwell 
1974) to solar photons in the near-UV is quite low, 
with a minimum ~0.25 near À % 2800 A. Unit optical 
depth for pure (conservative) Rayleigh scattering by 
H, for 2800 Å photons is reached at a level with 
P x 0.3 bar in the lower mesosphere. No simple 
molecules present near or above this level absorb in 
this wavelength region (NH is effective only for 
à < 2300 A), and complex molecules or very small 
absorbing dust grains have been invoked to provide 
the absorption (Axel 1972). Absorption in this region 
can be provided by a range of not implausible organic 
molecules (Sagan 1968; Khare and Sagan 1973). A 
mass absorption coefficient near 2800 A of ~(2Nore) 7? 
~ 108 cm? g~! would be required if organic molecules 
(beyond ethane) produced by Lea photons are re- 
sponsible. This is within the range of measured 
absorption coefficients for complex organic molecules 
produced by ultraviolet light in crude Jupiter simula- 
tion experiments (see, e.g., Khare and Sagan 1973). 


b) Cloud Models 


The total thermal infrared flux emerging from 
Jupiter gives Tert % 127 K. Most of this flux emerges 
at wavelengths A % 20 um to 60pm; the opacity 
source is mainly H, (with a small contribution from 
gaseous NHs); and most of the radiation emanates 
from levels with P ~ 0.7 bar to 0.2 bar (T ~ 150 K to 
~115K). Only a small fraction of the thermal flux 
emerges at wavelengths à < 15 um, but the opacities 
of the main gaseous constituents of the Jovian 
atmosphere are low at some of these shorter infrared 
wavelengths. In the absence of any particulate matter 
(or complex molecules), radiation with A ~ 8.2 um to 
9.5 um would emerge from levels with T ~ 155 K, and 
radiation with 4 ~ 5 um from even deeper levels with 
temperatures up to T % 300 K. 

Assuming solar abundances of all elements in the 
Jovian interior and assuming thermochemical equi- 
librium, one can calculate condensation temperatures 
Teon for various species (cf. Lewis 1969a; Weiden- 
schilling and Lewis 1973; Prinn and Owen 1976). Toon 
determines the level of each ‘“‘cloud bottom,” and 
above this level the gas-phase abundance is given by 
the saturation vapor pressure. The abundance and 
size distribution of droplets (or crystals) in the upper 
layers of a cloud cannot be calculated without detailed 
knowledge of the dynamics (see § III). The cloud 
bottoms for ice (and aqueous ammonia) are at levels 
with T ~ 280 K, for NH,SH at T ~ 200 K, and for 
NH; at T x 150K (P x 0.7 bar). 

Measurements of spectral reflectivities for solar 
radiation in the visible range are now available 
separately for the bright zones and for the darker 
belts in the Jovian equatorial and temperate latitudes 
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(Orton 1975c). The fraction of the sunlight that is 
absorbed varies from only ~0.2 in zones to ~0.3 in 
belts for red light, but varies from ~0.3 in zones to 
~0.5 in belts for blue light. The data at A ~ 20 um and 
45 um, separately for zones and belts (Orton 1975c), 
are compatible with a model which includes an 
optically opaque cloud top at the level (T x 145 K, 
P x 0.6 bar) for zones, but not for belts. This model 
is not unique, but the postulated cloud top occurs close 
to the level at which solid NH; crystals are expected to 
precipitate out if the zones represent the upward draft 
in the overall circulation. The spectral data for A x% 8 
to 14um (averaged over some zones and belts), 
together with spatially resolved data at A = 8.11 um 
and 8.45 um (Orton 1975b), indicate the presence of 
some obscuring matter at altitudes above the T = 
150 K level. The absorption spectrum of solid NH, 
crystals fits these data quite well, with more absorbing 
material over zones than over belts suggested by the 
data. Nevertheless, a small amount of absorption 
seems to be required over the belts as well, with an 
extent of only a few kilometers at levels with T % 
145 K to 150 K indicated for this thin haze. 

One question of primary concern is the nature and 
location of the red chromophores which are responsible 
for the absorption of visible sunlight over the belts (or 
at least for the preferential absorption in the blue). 
Upper limits to the location altitude can be inferred 
from the common observation that the red chromo- 
phores are sometimes overlain by time-variable white 
clouds. The relative heights are apparent because the 
interface between white and red material is frequently 
convex outward from the white clouds. As discussed, 
the highest lying clouds which are possible on Jupiter 
are ammonia cirrus clouds at levels with T ~ 145 to 
150 K. The red chromophores must therefore reside 
mainly at deeper levels with T > 145 K and P > 0.6 
bar. Emission at A % 5 um reveals infrared hotspots 
in only some regions of the belts (but in none of the 
zones) where the radiation emerges from deep levels 
with temperatures up to T % 300 K. The red colora- 
tion is anticorrelated with these infrared hotspots 
(Keay et al. 1973; Westphal, Matthews, and Terrile 
1974), which implies either that (a) the red chromo- 
phores reside above the level of some intermittent 
clouds which are opaque to 5 um radiation and which 
reside somewhere above the T ~ 300 K level; or (b) 
that the red chromophores are themselves intermittent 
and opaque at Sym. The infrared hotspots are 
associated with visually dark patches of bluish hue, 
probably connected with Rayleigh scattering rather 
than with blue chromophores, placing the blue regions 
at pressures of several bars, and T < 280 K (see Sagan 
1971). To summarize: The red chromophores reside 
mainly in the troposphere (T > 145 K), extending 
downward from the tropopause (T ~ 145 K) to at 
most the T % 300 K level. 

The condensation temperature of solid NH; is only 
slightly higher than the temperature (T ~ 145 K) 
expected for the tropopause, so that the solid NH, 
clouds reside in the upper layers of the troposphere. 
The pattern of patches of white clouds overlying the 
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red chromophores changes on a time scale of a few 
weeks, This is comparable to the total tropospheric 
eddy diffusion time ¢, in equation (3) if € ~ 1, but is 
also compatible with the longer total circulation time 
tor in Gierasch’s model, since time scales for patches 
of horizontal extent d, « D,, (the separation between 
zones and belts) are only of order (d,/Dz»)tor. The 
postulated intermittent 5 um clouds underlying the red 
chromophores could in principle reside anywhere 
between the T ~ 145 K and T ~ 300 K levels, but on 
Gierasch’s model these levels in belt regions represent 
a.continuous downdraft of warming material which 
could not form NH,SH or H,O crystals. Solid NH; 
clouds can be present in belts to varying degrees, 
because the upper horizontal circulation path from 
zones to belts can include layers just below the tropo- 
pause which carry NH, crystals. On this picture the 
red coloration would have to be provided by red 
chromophores which more or less coexist with NHs 
clouds (P ~ 0.5 to 0.7 bar), but with varying amounts 
of (i) high-level clouds with small NH, crystals, (ii) red 
chromophores, and (iii) clouds at lower levels with 
larger NH; crystals. 


c) Radiative Transfer Excursion 


Besides the red chromophores we are also interested 
in inferred near-UV chromophores: The measured 
albedo of Jupiter (Wallace, Prather, and Belton 1974; 
Savage and Caldwell 1974) to solar photons in the 
near-UV (A x 2000 Å to 3500 Å, say) is quite low, 
with a minimum albedo ~0.25 near A % 2800 A (after 
rising slowly toward shorter A, the albedo again seems 
to be lower for A < 2300 A). Absorption by HS, 
residing in the NH,SH clouds at lower levels has also 
been considered (Prinn 1970, 1974; Lewis and Prinn 
1970; Sagan and Khare 1971b; Lewis 1976). A 
computational problem is often presented by absorbers 
whose abundance decreases very rapidly with increas- 
ing altitude (because of condensation); we briefly 
review the relevant radiative transfer theory. 

Consider plane-parallel geometry and assume that 
the scattering is isotropic. Let 7 be the total extinction 
optical depth at some level (scattering plus absorp- 
tion), @ the single scattering albedo (ratio of scattering 
to extinction cross section) at the same level, and 
k = [3(1 — &,)]'. If « is independent of 7, the 
Eddington approximation is known even for slabs of 
finite thickness (Chandrasekhar 1960; Irvine 1975) and 
is discussed in our Appendix. The mean intensity J(7) 
as a function of depth is particularly simple for a 
semi-infinite slab of constant « (see eqs. [A3] to [A5] 
and eq. [A17]). Unfortunately, « can be a very steep 
function of 7 as in the example in Figure 2, which 
corresponds very roughly to absorption by HS and 
Rayleigh scattering by H of 2500 Å wavelength 
photons in the work by Prinn (1970). Prinn used 
equation (A3) for J(7) in the integrand of an integral 
which evaluates the absorbed energy layer by layer. 
Unfortunately, this use of equation (A3) can lead to 
gross errors when « is varying rapidly, and a more 
accurate numerical evaluation is given in our Appendix 
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Fic. 2.—A model radiative transfer calculation for an 
assumed form of the parameter «x, as a function of total 
optical depth 7. Numerical results for the mean intensity J are 
given in the solid curve. The dashed curve is an analytic 
approximation, discussed in the text. 


for the example in Figure 2. Fortunately, the correct 
result for the fraction f of the incident energy which is 
absorbed is fairly close to the following simple 
approximation, fọ: Let To be the total optical depth at 
the level where «x equals some predetermined number, 
say x = 0.5 (in our example 7, % 7.5). Let fo be the 
fractional energy which would be absorbed according 
to the Eddington approximation if one had no 
absorption at all for 7 < 7) and a perfect absorber at 
Tọ For diffuse incident illumination (instead of 
averaging over different angles of incidence) one has 
So = (0.7579 + 1)~1 if To > 1 (see eq. [A16]). For the 
numerical example in the Appendix the value of fo 
corresponding to the choice «x = 0.5 for ro fits the 
correct f quite well. 


d) Near-Ultraviolet 


Solar photons in the near-UV (A = 2000A to 
3000 A, say) are much more abundant than Le 
photons. It is of interest to speculate about materials 
which absorb in the near-UV to see (i) if their absorp- 
tion can explain the low observed albedo in the near- 
UV, and/or (ii) if the end products of reactions 
initiated by such absorption can provide the red 
chromophores for the troposphere, and (iii) how much 
complex organic material can be so produced. 

The absorption of solar radiation with à ~ 2300 A 
to 2700 A by H3S in the upper layers of the troposphere 
has been discussed by the authors listed above (for 
À > 2700 A the absorption by HS is too weak and 
photons with A < 2300A are already absorbed by 
NH; at higher altitudes). The total flux in this wave- 
length range (diurnal mean) is about 1 x 10** photons 
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cm~?s~?, and some fraction f, of these photons are 
absorbed by HS. We first present estimates if an ad 
hoc value of f, = 0.1 is assumed (and £ = 1 in eq. [3]). 

Laboratory simulation experiments by near-ultra- 
violet irradiation of CH,, C,Hs, NH3, H2O, and HS 
mixtures at ~1 bar pressures have been performed, 
but without initial excess hydrogen (Sagan 1971; 
Sagan and Khare 197la, b; Khare and Sagan 1973, 
1975). Hydrogen sulfide is the primary photon acceptor 
in these experiments. Amino acids were produced in 
these experiments with a quantum yield of ~10~* for 
molecular weights ~ 100 (~ 10-2 amu photon ~*+). The 
quantum yield of all organics in such experiments is of 
the order of 100 times larger (of the order of 1 amu 
photon~*). Some red chromophores were abundantly 
produced, consisting mainly of polymeric sulfur, but 
including a substantial fraction of organic compounds. 
The overall mass absorption coefficient for blue light 
of these chromophores is k ~ 10° cm? g~?. If such 
particles are produced in Jupiter at the rate of one 
S-bond per absorbed photon (R ~ 1072S cm~2s71 ~ 
5 x 107" gcm~?s~?), their column density, N, over 
a few scale heights down to the thick water clouds, is 
needed. According to equations (3) and (8), N z 
2RH?/K, ~ 107+ gcm~? and the optical depth to 
absorption in the blue is 7 = kN ~ 0.03. This is 
almost two orders of magnitude smaller than the 
observed optical depth. Moreover, the reaction of 
eight sulfhydryl radicals to form Ss is poisoned by 
competing reactions, and the Sg quantum yield must be 
significantly less than 1. 

In the laboratory experiments cited, organic mole- 
cules, but not polymeric sulfur, were produced from 
superthermal hydrogen atoms. In the presence of a 
thousand times more H, (by number), the quantum 
yield might be expected to decrease by about the same 
factor, if energy loss of the fast atoms is mainly by 
elastic collisions, to ~10~* amu photon~?. The total 
production rate would then be quite small, R ~ 5 x 
10714 gcm~?s~? (which is still 10 times larger than 
our estimated R,,, for the mesosphere). The corre- 
sponding optical depth for organic chromophores is, 
then, many orders of magnitude below the observed 
value. Many uncertainties remain, however, including 
indirect effects of scattering. In the laboratory ex- 
periments cited, large quantities of molecular hydrogen 
are produced during the irradiation because the end 
products are unsaturated, and it may be that the 
correction factor from laboratory to Jovian situation 
is closer to 107? than to 107°. It is also possible that 
somewhat longer wavelength ultraviolet light—perhaps 
as long as 2900 A—is still effective in these syntheses. 
Recent experiments (Ferris and Chen 1975) have 
shown that 1849 A irradiation of mixtures of CH, and 
NH, in a 10:1 excess of Hz produce a variety of 
organic molecules in very high yields, presumably by 
the generation of hot hydrogen atoms from ammonia 
photodissociation. These results suggest that the effects 
of H, dilution on poisoning organic photochemistry 
may be slower than linear. Hence the contribution of 
organic polymers to near-UV photoproduced optical 
frequency absorbers on Jupiter may be significant. 
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Polymeric sulfur is not produced by hot hydrogen 
atoms in such experiments, and so the scaling from 
laboratory to Jupiter is more direct in this case. How- 
ever, even here, and even under the excessively optimis- 
tic assumption that every H,S photodissociation event 
leads to polymeric S, the implied optical depth falls 
short by two orders of magnitude from matching the 
observed values; moreover, pure polymeric sulfur fits 
the observed optical properties of the Jovian red 
chromophores poorly at best (Rages and Sagan, 
in preparation). 

We have seen that a slow, deep circulation pattern 
is likely to extend all the way from the tropopause 
down to the level with T ~ 300 K (Gierasch 1976) 
with £ ~ 0.01 in equation (3). Abundances of reaction 
products are proportional to €~/ and also to ({,/0.1), 
which was claimed to be as large as 4 by Prinn (1970). 
One might then expect HS photolysis to dominate 
chromophore production, but this still seems unlikely 
for three reasons: (1) As discussed in § [Vc and in our 
Appendix, Prinn’s (1970) estimate of f % 0.4 is likely 
to be an overestimate since fo = (0.757) + 1)? is 
smaller. Thus f, < 0.2 is more likely even in the ab- 
sence of any scatterers (other than H) or other 
absorbers above the H.S levels. The smallness of f, is 
also suggested by the fact that the observed albedo of 
Jupiter does not increase from A ~ 2800 A (where 
HS is ineffective) to À ~ 2600 A (where HS absorp- 
tion should already be strong). (2) Since HS cannot 
itself be the cause for the low albedo in this wave- 
length range, whatever absorbers are the cause must 
lower f, further. (3) Most important, the systematic 
deep circulation pattern mentioned before depresses 
H.S photolysis both in zones and belts: In the updraft, 
represented by zones, HS is abundant in the NH,SH 
clouds, but NH, clouds are above these layers and 
prevent most of the solar radiation from penetrating 
down. In the downdraft, represented by belts, the 
NH, clouds are absent (or at least diminished), but 
the material moving downward is practically devoid of 
H,O and H.S (which condensed out on the previous 
upward journey). HO and HS are again abundant in 
deep levels with T > 300 K, but little of the near-UV 
penetrates to such depths. 

NH, absorbs solar photons of wavelength A < 
2300 A, which are somewhat less abundant than 
photons between 2300 A and 2700 A. On the other 
hand, NH, should not suffer the shielding difficulties 
experienced by H.S because it resides in the highest 
cloud layer which extends up to the tropopause. NH3 
should therefore be abundant in both zones and belts 
and even flows into the lower mesosphere at a non- 
negligible rate. It is therefore likely that most of the 
photons with A ~ 2300 A are absorbed by ammonia 
(Strobel 1973b), which is compatible with the observed. 
decline in albedo from A > 2300 Å to A < 2300A. 
Hydrazine particles are also produced by ammonia 
photolysis and may provide a suitable explanation for 
the low albedo near A ~ 2800 A (Prinn 1974). More- 
over, organic molecules produced by NH, photo- 
dissociation at A < 2300 A could conceivably be the 
principal source of red chromophores, if the quantum 
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yields of such reactions with 108:1 H, dilutions are as 
high as the experiments of Ferris and Chen (1975) 
suggest. 


e) Optical Frequency Photoproduction of Chromophores 


We here consider two possibilities: If the circula- 
tion in the upper part of the troposphere is as slow as 
Gierasch (1976) suggests, then € < 0.01 in equation 
(3) and column densities obtained from equation.(8) 
are increased by two orders of magnitude; thus 
polymeric sulfur or complex organics might provide 
the red chromophores. But, as we have seen, there are 
other problems with the contention that the chromo- 
phores are produced by H2S photolysis. These problems 
are avoided if shorter wavelength photolysis of NH, 
is an adequate source of complex organic chromo- 
phores, but this is a questionable proposition consider- 
ing the lower photon flux at A < 2300 A, and whatever 
is the correct value of the poisoning of organic photo- 
chemical synthesis in a 108:1 H, dilution. If, on the 
other hand, é turns out to be of order unity, then 
ultraviolet photochemistry cannot in any case provide 
the observed chromophores at a sufficient rate. Thus 
in both cases we are led to the possibility that photo- 
chemistry at optical frequencies, where more total 
energy is available, is involved in chromophore pro- 
duction. But no chemical bonds of even moderately 
abundant molecules can be broken by photons 
of wavelengths longer than 3500 A on Jupiter, from 
which it follows that optical frequency photochemistry 
must be at least a two-photon process. However, 
such processes appear to be exclusively biological. 
We therefore find ourselves led unexpectedly to the 
hypothesis that the Jovian chromophores are biological 
in origin and that there is an abundant biota in the 
Jovian clouds. 

Let us assume that the optical chromophores are 
contained primarily within Jovian organisms in the 
upper troposphere—organisms driven by visible light 
photosynthesis, perhaps primarily in the blue, and 
utilizing the abundant methane, water, and ammonia. 
Terrestrial algae (using CO, instead of CH,) under 
laboratory conditions (Lehninger 1972) have a quan- 
tum yield for synthesizing glucose of about 4 amu 
photon™? (in the visible). The equivalent maximum 
production rate of biological material in the Jovian 
troposphere would be about R = 3 x 10-° g cm~? 
s~1. With a mass absorption coefficient of k ~ 10° 
cm? g~t, such organisms would yield an optical depth 
of 7 ~ 2, slightly more than required by the observa- 
tions. If the parameter é turns out to be of order 10~?, 
and the red chromophores are due to Jovian photo- 
synthetic organisms, R and/or K would have to be 
smaller. 

Thus Jovian organisms having metabolic and photo- 
synthetic parameters typical of terrestrial algae are 
able to account for the optical chromophores on 
Jupiter. Such organisms must, however, be adapted to 
the Jovian environment, and in particular must 
maintain a steady-state population in the face of 
convective pyrolysis. Likely ecological niches for such 
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organisms, as determined by the hydrodynamics on 
Jupiter, are discussed in the following sections. The 
following discussion can, however, be treated inde- 
pendently of the hypothesis that Jovian chromophores 
are biogenic. 


V. SINKERS AND FLOATERS 


If there is a Jovian biology, we would expect it to 
fill a rich variety of ecological niches. The best 
terrestrial analogy seems to be the surface of the sea. 
Oceanic phytoplankton inhabit a euphotic zone near 
the ocean surface where photosynthesis is possible. 
They are slightly denser than seawater and passively 
sink out of the euphotic zone and die. But such 
organisms reproduce as they sink, return some 
daughter cells to the euphotic zone through turbulent 
mixing, and in this way maintain a steady-state 
population. (The early stages of sinking move the 
phytoplankton from a region of depleted resources 
into a region of abundant nutrients, thus stimulating 
replication.) 

A more elaborate adaptation is provided by fish and 
other organisms with float bladders which use meta- 
bolic energy to maintain a habitat at suitable pressure 
levels. They are generally not photosynthetic auto- 
trophs, but rather heterotrophs living on organic 
molecules produced by autotrophs. A third ecological 
niche is filled by marine predators, one step further up 
the food chain, which hunt heterotrophs. 

In the following discussion we will consider three 
comparable ecological niches on Jupiter: The primary 
photosynthetic autotrophs, which must replicate 
before they are pyrolyzed, will be described as sinkers. 
A second category of larger organism, which may be 
either autotrophs or heterotrophs but which actively 
maintain their pressure level, will be described as 
floaters. A third category of organisms actively seek 
out other organisms; we call these hunters, although, 
as we shall see, the distinction between hunting and 
mating under these conditions is not sharp. Finally, 
there is a category of organisms which live almost at 
pyrolytic depth. They are scavengers, metabolizing the 
products of thermal degradation of other organisms. 
For our purposes, these pyrolytic scavengers are 
identical with floaters. In the following two sections we 
estimate the growth time ¢,, required for an organism 
to double its mass. 

We first consider passive organisms with positive 
excess density Ap, above the gas density p,. Such 
organisms falling under gravity through the atmosphere 
rapidly acquire a vertical drift velocity, wz, which 
depends on the location and size of the organism and 
increases with increasing Ap,. There is a premium on 
keeping Ap, and hence w; small, and we consider only 
organisms in the shape of thin, gas-filled balloons. 
Such organisms have been briefly proposed before, 
both for Jupiter (Shklovskii and Sagan 1966) and for 
Venus (Morowitz and Sagan 1967). 

For simplicity we treat a spherical shell of outer 
radius a and skin thickness d « a. In general we con- 
sider organisms filled with gas of the same temperature 
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and mean molecular weight (2.2) as the surrounding 
ambient atmosphere. The lack of buoyancy is then 
controlled by the excess averaged density due to the 
skin, Ap, = (3d/a)p,, where we assume the biological 
skin material to have unit density (1 g cm~). 

For sufficiently small radius a of an organism, the 
Stokes formula in equation (11) holds (with par —> 
3d/a) for the drift velocity w4 and the drag coefficient 
Cp & 12v/awg is large. Let are be a critical radius for 
which this relation gives Cp = 1 (Reynolds number = 
24), so that 


_ 2 gad 
Ya = 3 vp, 


are = (6r°p,/gd)™" . (17) 


For large organisms with a > ag, (large Reynolds 
number), the flow around the body becomes turbulent 
and the drag coefficient Cp (= drag force/0.57a?p,w”) 
is Close to 0.5 (see, e.g., Batchelor 1970). In this regime 


wa © A(ed]p,)!? . (18) 


We shall have to consider separately organisms 
living below and above the tropopause. The tropo- 
sphere provides conditions closest to terrestrial ones 
(similar temperatures and abundant water in liquid or 
vapor form). We assume that the organisms can 
survive there over a range L of vertical height (L ~ 5H, 
say), centered approximately on the level with T ~ 
300 K. An important time scale is then a nominal 
drift time t4 = L/wg. We consider the skin thickness d 
as an unknown parameter, with d ~ 1074 cm possibly 
a typical value (biological unit membranes on Earth 
have d < 10~° cm, but mechanical strength is required 
for the skin). Using equations (1) and (11) and 
A a coefficients at the level with T = 300 K, we 

nd 


for a K are, 


L ( T ys 1 cm 10-* cm 


w — [— 4 
ta X oH 300E a 7 1.3 x 10*s 
for 


-4 
G2, © (2 7m sek V 0.039 com. (9 


For a > are, on the other hand, equation (18) gives a 
drift time which depends on d but not directly on 


radius a, 
-4 
al A yr mas x 105s. (20) 


la X 3 (300K d 


The eddy diffusion coefficient K, in the convective 
troposphere is fairly large, and we also have to con- 
sider the contribution of eddy diffusion to the organ- 
ism’s downward motion over a vertical distance L. We 
are interested in cases where L is larger than H and 
the effective circulation time ty, for downward 
diffusion is then not given by L?/K,: With L > H and 
K: a slowly varying function of height, the first 
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derivative term in equation (7) can be neglected and 
eddy diffusion has the same effect as increasing w, to 
(wa + H~1K). The relevant time scale is then linear 
in L and is approximated by the smaller of t4 and the 
quantity 


LH = (300 


à 300K 


sips 3.7 x 10° 
for = Fle = K, ~” 5A 1X S, 


(21) 


where we have used equation (3) and assumed é ~ 1. 
If d x 10-* cm, then eddy diffusion competes with 
gravitational fall only for microorganisms with 
a = 30 um. Thus Jovian organisms the size of small 
terrestrial protozoa and prokaryotes have typical 
times for falling through the troposphere to pyrolytic 
depths of 1 to 2 months. To maintain a steady-state 
population, they must only replicate in that time scale. 

In the Jovian mesosphere the temperature is low (we 
assume T ~ 150K) and water is absent, but some 
biologies have nevertheless been envisioned for such 
environments (see, e.g., Pimentel et al. 1966; Sagan 
1970). We shall see that the lower, denser layers of the 
mesosphere are the most advantageous for the growth 
of such organisms; eddy diffusion is negligibly slow 
here compared with the downward drift velocity wa. 
For organisms operating between a lower pressure P 
and the tropopause (pressure ~1 bar), the drift-time 
L/w, becomes 


x 930s 


neta 1 bar\ 1 cm 10-4 cm 
iad P a d 


for 


-4 
P (2 a l r) “0.044 cm. (22) 


For larger organisms, a > dpe, the expression for 
large Reynolds number gives 


EA 1 bar\(10-*cm P 
= P d lbar 


To obtain the pyrolysis time scale this number must be 
added to the appropriate tropospheric time scale given 
by equations (19), (20), and (21). 

The expressions given above for tz hold only for 
balloon organisms with ambient atmosphere inside. 
For organisms capable of pumping gas, buoyancy can 
be achieved while maintaining pressure equilibrium by 
keeping the interior gas pure hydrogen with molecular 
weight u = 2.0 instead of the ambient hydrogen- 
helium mixture with p = 2.2. A pumped organism can 
thus float if Ap, = (3d/a)p, < 0.1p,. This requires 
radii a larger than a threshold radius a, for floating, 
given by 


d T 1 bar 
an X (tea) (55x) (=) a doem UA 


1/2 
) 1.6 x 10*s. (23) 
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Instead of maintaining a hydrogen interior, an organ- 
ism could use a fraction of its metabolic energy release 
to heat its interior gas. The temperature difference AT 
between interior and exterior is determined by the 
metabolic rate and the thermal diffusion coefficient of 
the gas. Except at very high pressurés, the thermal 
diffusion is fast enough to keep A7/T « 0.1, so that 
smaller values of Ap, and larger radii than a, are 
required. Buoyancy from a hot interior is thus likely 
to be less advantageous than pumping out helium. 
Note that Jovian floaters are macroscopic organisms 
(eq. [24]). Specialized floaters, with a range of organ 
systems, might have an effective value of d ~ 1 cm as 
terrestrial reptiles and mammals do; in that case, 
floaters would have kilometer dimensions, a size 
within the resolution capability of the Mariner 
Jupiter/Saturn flyby imaging system. The existence 
of very large floaters is, however, limited by metabolic 
constraints described in the next section. 


VI. GROWTH OF PASSIVE SINKERS AND FLOATERS 


We defer for the moment further questions of 
buoyancy and of powered locomotion, and consider a 
range of metabolic niches. Organisms—sinkers or 
floaters—may be photosynthetic autotrophs, obtaining 
their free energy directly from sunlight. They may also 
be heterotrophs, of which we distinguish two classes: 
(i) those which feed on organic molecules of non- 
biological origin (e.g., ethane, an ultraviolet photo- 
product), which reach the organism by molecular 
diffusion; and (ii) those which feed on organic matter 
in smaller organisms or their fragments which reach 
the larger organism because of the difference in 
vertical drift velocity wą between predator and prey. 
Active hunters described in the following section are a 
related subcase. In all cases we must estimate the 
growth time ¢,, = [2d ln a/dt]~* and compare it with 
the drift time ¢, (or with tar if that should be smaller). 
We shall see that the ratio ¢,,/tg increases with the 
radius a of the organism; we are interested in the 
maximum radius a,,,, for which this ratio is unity. If 
max is appreciably larger than the minimum size amin 
for a particular life style, then a biology can be 
maintained by an organism fragmenting into many 
smaller organisms or dispersules before it reaches size 
max, 1.€., before it has drifted downward to a pyrolytic 
death. If ta > ty, for the dispersules, they can be 
circulated upward, grow to a,x, and complete the 
life cycle. 

We saw in §IV that the most favorable region for 
direct photosynthesis is probably the upper tropo- 
sphere, just above the water clouds (with T ~ 300 K, 
P x 5 bar). Methane, ammonia, and water building 
blocks are abundantly available there to a photo- 
synthesizing organism, and the growth rate is con- 
trolled by the deposition of free energy from solar 
photons. Furthermore, the observed optical chromo- 
phores can be explained by such organisms at such a 
level. We saw that a value for the production rate of 
R x 3 x 107° g cm~? s~! seems reasonable, and we 
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find for the growth time 


(25) 


Equating this time with ¢, in equation (19) then gives 
for dmax, the maximum radius to which an- organism 
can grow, 


Amax © (10-74 cm/d)?0.1 cm. (26) 


Values of amin $ 107? seem perfectly feasible, so that 
one organism can produce more than 10* dispersules 
for which t, «4x 106s % tar The steady-state 
population of such organisms would then be con- 
trolled only by competition for sunlight, and not by 
gravitational fallout to pyrolytic depths. 

Consider next the diffusion onto the skin of an 
organism of some organic molecule with molecular 
weight u and abundance in the surrounding atmosphere 
of po (in g cm~? of atmosphere). Ethane is the lowest 
mass (u = 30) organic molecule which carries free 
energy (relative to methane, the prevalent carbon 
molecule at thermal equilibrium). The diffusion co- 
efficient D of ethane in the Jovian H-He mixture is 
about 0.4 times the viscosity v of this mixture (Strobel 
1973), and for heavier molecules D scales reciprocally | 
as the cross section of the molecule. For approxi- 
mately spherical organic molecules we adopt 


D x 0.4r(u/30)-2!3,, (27) 


with v given by equation (11). If v > D > (aw,/12), 
the diffusion is almost independent of the | motion of 
the sphere, and the steady-state rate for diffusion onto 
the surface is 47a?(p)D/a). Let e be the ratio of in- 
creased body mass of the organism to the mass of 
ingested metabolites. Assuming unit density for the 
skin material, the growth time is then 


ter = ad]2epoD . (28) 


At a mesospheric level with pressure P (in bars), for 
molecules produced at a mass rate of R = pO at 
some higher level, equations (6) and (9) give po ~ 
(u®H/100)P In P+. The growth time in the meso- 
sphere is then 


a d (u./30)?° 
‘ex © Tom 107? cm In (1 bar/P) 


-14 
x (>) 0.03 18 x 108s. (29) 


For ethane an efficiency factor e ~ 0.03 seems reason- 
able and R ~ 5 x 10714 gcm~?s~? has already been 
discussed; for UV photoproduction of more complex 
organic molecules we may leave e ~ 0.3 and R ~ 5 x 
10-1 gm cm~? s7}, so that Re is the same, but the 
factor u” makes heavier molecules less favorable. 
Replacing the slowly varying factor In (1/P) by 3 and 
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Re by 1.5 x 10715, we find from equations (22) and 
(29), 


Amax ~ (1074 cm/d)(u/30)-" x 0.07cm. (30) 


If Me is particularly large and/or d particularly 
small, an organism can grow beyond size are, and the 
diffusion rate now depends on the drift speed w4, given 
by equation (18). We omit the complicated transition 
region a X% dp, and make only order-of-magnitude 
estimates: There is a boundary layer in the flow 
around the sphere of thickness b © a(4v/aw,)'!? « a, 
such that the flow speed a normal distance y from the 
body is ~(y/2b)w,g. With D < v the typical stand-off 
distance for diffusion onto the body surface is bp ~ 
2b(D/v)"8 and the diffusion rate is ~ 47a7(p)D/bp), so 
that 


ter ~ (2d/epo)(a/ Dwa)? 0| DY" . (31) 


For the mesosphere this gives 


a \12 d 3/4/1 bar\ 1/4 (u/30)*/? 
ter ~ (15) (sea) (3 in /P) 


5 x 10-14 0,03 , 
x (SE B s x 10 S. (32) 


The same basic formulae apply in the troposphere, 
but pọ is now much smaller for the same value of ® 
because of the larger eddy diffusion coefficient, K,. On 
the other hand, we have seen that photosynthesis may 
be going on there, and this may release complex 
organic molecules at a faster rate ®. Substituting 
po = RH/K, and using equation (3), we find (for 
a < dre) 


~a d T p\23 
sr ~ Í cm 10-4 cm \300 K 30 
3 x 107° 0.3 i 
x ( D 3)8 x 10 s. (33) 
Comparison with equation (19) shows max to be of 
order (10-74 cm/d)(30/4)*80.013 cm, which is just 
macroscopic. 

If small tropospheric photosynthetic organisms have 
radii a, < 0.01 cm, their drift speeds are even smaller 
than the eddy circulation speed and their abundance 
is still given by pọ = RHA/K,. A larger organism of 
radius a, drifting downward with speed w4, can grow 
by coalescing with such small organisms if the en- 
counter satisfies certain conditions. These conditions 
are similar to those governing the coalescence of rain- 
drops (Mason 1972; Rossow 1976). In analogy with 
equation (13) one finds, independent of a,, and pro- 
vided d «a< a, the requirement that a> 
(10- * cm/d)(7/300 K)0.002 cm. If these conditions are 
satisfied, the larger organisms incorporate mass from 
the smaller ones at a rate of 4ra?wapo, so that ter X 
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(4d/R)(K,/Hw,). If a > are ~ 0.04 cm, then 


a d Wa T \ (3 x 10-° 0.3 
i 10-* cm 300 ® € 


x 2.5 x 10's. (34) 


Neither ¢,, nor tz in equation (20) depends explicitly on 
radius a. An organism can then grow indefinitely as 
long as (d/10~* cm) < (7/300 K)*'?. 

Thus, provided the membrane thickness stays 
smaller than about 1 um, such passive sinkers can 
grow by coalescence with smaller organisms to sizes 
larger than an as given by equation (24). This provides 
one evolutionary pathway from sinkers to floaters. 
However, the condition on membrane thickness is 
very restrictive. Floaters of this sort are unlikely to be 
much more stable than soap bubbles, and the organ 
systems which large organisms require could not be 
accommodated with the thin skin thickness. However, 
as we shall see in the following section, the hunting 
lifestyle provides a way out of this dilemma. 

To summarize the situation so far: Sinkers can grow 
and produce reproductive dispersules in a stable, 
steady-state life cycle if they either (a) are mesospheric 
heterotrophs or photoautotrophs, or (b) are tropo- 
spheric photoautotrophs. Passive floaters exist only if 
(1) their skin can be made particularly thin, and (2) a 
significant fraction of the smaller photoautotrophs are 
available as food. 


VII. POWERED LOCOMOTION AND COALESCENCE 


We saw that the efficiency of utilization of meta- 
bolites diffusing onto the surface of a freely falling 
organism decreases with increasing mass of the 
molecule. However, for even more massive metabolites 
and organisms approaching each other with some 
relative drift velocity, growth by coalescence is again 
moderately efficient. The growth rate by coalescence 
can be increased greatly if organisms are capable of 
powered locomotion and can steer toward each other. 
We shall not discuss the predator-prey relationships of 
a strict hunting doctrine, but rather only an ecology in 
which an organism approaches another, coalesces into 
one larger organism which approaches other organ- 
isms, etc. It is clear that this hierarchical mating 
doctrine can also serve the traditional terrestrial 
function of mating, namely, the exchange of genetic 
material. A strict hunting doctrine in which one of the 
partners in such a coalescence is noncooperative will 
be less efficient than the doctrine discussed here. For 
convenience, however, we will describe both doctrines 
as hunting. In the following discussion we will assume 
that only one partner in such a mating event is active 
in the search; for this reason mating, hunting prey, and 
hunting abiological organic molecules are not ex- 
tremely different. The doctrine in which both organ- 
isms are active in seeking partners would be more 
efficient than the ones described below. 

For simplicity we assume the presence, with number 
density n, of organisms with radii of order a, at some 
upper injection level z = Z,. The organisms coalesce 
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and grow as they drift downward, so that the size 
distribution changes with height z. Let n(a, z) and 
p(a, z) = (47a*d)n(a, z) be their steady-state number 
density and mass density, respectively, per logarithmic 
interval of radius a. In principle we have a relation like 
equation (7) separately for each value of a, with the 
right-hand side representing the flux from one size 
range to another. We shall make only order-of- 
magnitude estimates for ā(z), the mean radius of 
organisms at level z, and for n(z) = f d In an(a, z). 

We assume that the organisms are sufficiently large 
that w,(z), the drift velocity for radius a(z) at height z, 
is much larger than K/H in equation (7). Multiplying 
by mass (assuming as usual unit density for the skin 
material) and integrating over sizes, we find 


[47ra~#(z)d]w,(z)n(z) = €(z)R. (35) 


In this relation R is the constant rate in g cm~? s~! at 
which mass in the form of organisms enters the upper 
level z = Z,; e(z) is a slowly varying factor (slightly 
smaller than unity) which takes account of any in- 
efficiency in the utilization of biological material due 
to waste or expenditure of free energy. We shall have 
to evaluate the growth time ¢,,(z) = (2d In a/dt)~* for 
organisms with a ~ ā(z) at level z. The function A(z) is 
then given in terms of ter by 


din aldz = [2t,e,w,(z)]-?. (36) 


For the integrated column density N(a) = f dzn(a, z) 
we also find N(a) % 2t,wgn(Z), where Z is such that 
a(Z) = a. The total mass contained in organisms in 
different size ranges is then given by 


(4ra*d)N(a) X 2tey(a)e(z)R . (37) 


Assume that an organism, while drifting downward 
with speed w4, can expend enough mechanical energy 
to acquire a horizontal velocity component of up to 
Owa, where 0 is some angle less than | radian. Its 
vertical speed relative to other similar organisms is 
about 0.5w,z and in time ¢ a volume (67/6)(w,t)? is 
accessible to a searching organism. When n(z) times 
this volume is unity, the organism has a good proba- 
bility of coalescing with a mate and doubling in mass. 
With n(z) given by equation (35), the growth time is 
then 


tg 


4 ( 3rã?d y. (38) 


j 5 h2 eRO?Ww:? 


If @ > age we have, independent of d both for the 
troposphere and the mesosphere, 


E? L a \2 175K P 3 x 10°°0.3 H3O s 
sr ~ | 62 \1 cm T 1bar R € ° 


(39) 


For such organisms the troposphere is quite favorable 
even in the absence of photoautotrophs as prey 
(R ~ 5 x 10714 gcm~*s~?); for instance, with T x 
300 K and d ~ 10-*cm, we have an ~ 6cm and 
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TABLE 2 
MAXIMUM HUNTER SIZES 


SOURCE OF METABOLITES 


Abiological 
0 Organic Matter §Photoautotrophs 
2 x 10°? radians..... 6cm 40m 
Pick beeen EEN 3m 2km 
(ere eee eee renee 30m 20 km 


Amax ~ (6/0.002)6 cm, so that buoyancy can be 
achieved with quite small values of 0. For R ~ 3 x 
10-°gcm~*s~1, ay ~ 6cm requires only amaz ~ 
(6/3 x 107°)6 cm; microradian maneuverability suf- 
fices. The maximum sizes of organisms embracing the 
hunter doctrine as a function of the maneuverability, 
6, and the source of food are given in Table 2. We see 
that very large hunters are permitted, within the 
resolution limits of the Mariner Jupiter/Saturn imaging 
system. 

To derive an explicit relation between typical radius 
a and vertical level Z(a), we must integrate equation 
(36) with ter given by equation (39). The result depends 
on whether ô is constant or a function of a. Neglecting 
the variations of 0, wg, and e gives (for a much larger 
than the original size a, at level Z,) 


Z, — Z(a) = 3wa(Z)ter(@) . (40) 


We still have to calculate the mechanical energy 
expended by an organism in maintaining its horizontal 
hunting speed Ow,. In pure downward drift the rate of 
energy dissipation equals the rate of gravitational 
energy release mgwa. The equivalent rate for maintain- 
ing the horizontal speed (with 0 < 1) is at most Omgwg. 
(With judicious sailing and soaring it could be appreci- 
ably less.) The relevant quality factor is the ratio 
AE/Am of mechanical energy expended in finding a 
partner to the increase in mass of biological material. 
(Biochemical energy storage is of order 10** ergs g~t 
and a mechanical quality factor ~10° ergs g7? is a 
very modest assumption.) The heat of combustion 
(full oxidation) of molecular hydrogen is about 10 
times that of glucose and more than 3 times that of 
such fats as stearic acid. It is not the unavailability 
of reduced compounds but rather the unavailability of 
oxidized compounds on Jupiter which probably sets 
limits on the efficiency of internal metabolism of 


hypothetical Jovian organisms. We have 
AE 
Am <= OgWater > (41) 


which depends quite weakly on temperature and 
pressure. For the upper troposphere 


AE (a \"(  d_\2 
Am ~ (; =) (isan) 

( 3 x 107°0.3 

x (0 ee 


1/3 erg 
6 
2 5 x 10°-* (42) 
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represents quite a modest power expenditure, except 
for the largest organisms. A 100m radius hunter of 
photoautotrophs with 0 = 0.01 radian excursion 
capability and a 1 cm skin thickness requires a quality 
factor of 107 ergs g-1. The same quality factor is 
required by a 1 m radius hunter of abiological organic 
matter with 6 = 0.03 radians and d = 1 cm. 

Even smaller power expenditures will be required if, 
as would be likely, evolution selects organisms with 
sensory systems which increase the efficiency of hunt- 
ing food—such as, e.g., optical sensors for detecting 
chromophores. Likewise, acceleration sensors capable 
of detecting convective ‘‘thermals”?” would ease the 
restrictions on floater sizes and energy expenditures; 
and habitats near the boundaries between upward and 
downward convective regimes would improve the cost 
of Jovian biology. 

The path from sinkers to hunters to floaters is an 
evolutionary path. It provides a possible sequence for 
the evolution of floaters which, as we have seen, 
would be difficult to understand directly from sinkers. 
Once floaters are established, they may be able to 
replicate without passing through the sinker stage at 
any point in their life cycles. For this to occur, of 
course, the minimum radius of their dispersules must 
be larger than an. 

While we have distinguished passive sinkers, 
hunters, and floaters in this discussion, it is clear that 
a single highly evolved Jovian organism might partake 
of qualities of all of these life styles—for example, at 
different stages in the life cycle. Indeed, it may be that 
only after the hunting doctrine is adopted by sinkers 
can growth to sizes large enough for the evolution of 
floaters become possible. We have in this discussion 
made no distinction among various locales on Jupiter; 
but it is clear that some locales—the Great Red Spot, 
for example—may be more favored than others be- 
cause of higher abundances of organic molecules, 
prevailing updrafts, or other reasons. 

Among other objects in the outer solar system which 
possess atmospheres, only Titan and Saturn exhibit 
colorations similar to those on Jupiter. Because of the 
absence of convective pyrolysis on Titan, the chromo- 
phores there can quite readily be ultraviolet photo- 
produced organic molecules (Sagan 1973). However, 
the physical environment of Saturn is very similar to 
that of Jupiter in the respects relevant to the arguments 
of this paper, and we therefore tentatively consider an 
airborne biota on Saturn as well. 

The test of these ideas is, of course, observational. 
Flyby, orbital, and entry probe spacecraft each can 
perform significant tests of these ideas. For example, 
entry gas chromatograph/mass spectrometers for the 
near future are anticipated to have detectivities of 
organic matter as good as about 10°7° gcm~3. The 
predicted steady-state density of complex organic 
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TABLE 3 
STEADY-STATE DENSITY OF ORGANIC MOLECULES (g cm 7°) 
Abiological Biological 
Lower Mesosphere........... 2 x 1071 9 x 1078 
Upper Troposphere.......... 2 x 10775 9 x 10-710 


molecules on Jupiter is RH/K, and is roughly estimated, 
using the numerical values of this paper, both for the 
lower mesosphere and upper troposphere and for 
abiological and biological sources of organic matter, 
in Table 3. While most biological organic matter will 
be in organisms and not in free molecules, the densities 
predicted in Table 3 are sufficiently large that a high- 
sensitivity entry gas chromatograph may bea significant 
test of the biological hypotheses of the present paper. 
Because of the slow convective velocities of the 
mesosphere, such instruments should be capable of 
working at such pressure levels. For comparison, the 
density of dissolved organic matter at 2000 m depth in 
the north-central Pacific is equal to 5 x 1077 gcm~3, 
and of particulate organic matter, 3 x 107? g cm~? 
(Williams and Carlucci 1976). 


VII. FINAL REMARKS 


The possible existence of indigenous Jovian organ- 
isms is also relevant to the question of sterilization of 
spacecraft intended for entry into the atmosphere of 
Jupiter. Even if the ambient environment is incon- 
sistent with the growth of microbial contaminants 
from Earth, the internal environment of Jovian 
organisms may be much more hospitable. The 
replication of terrestrial contaminants in the Jovian 
clouds also depends on the availability of trace 
elements. We have seen (§ III) that indigenous Na is 
probably missing from the upper troposphere, and 
that indigenous Mg must certainly be missing. 
Magnesium ions are essential] for nucleic acid replica- 
tion and a wide range of other biological functions; if 
Mg is missing from the clouds, the likelihood of 
terrestrial biological contamination of Jupiter becomes 
nil. However, exogenous sources of sodium, magnes- 
ium, and other trace elements, particularly from 
micrometeorites, while small, cannot be neglected 
(§ III). Consequently, it seems judicious not to exclude 
prematurely the possibility of biological contamina- 
tion of Jupiter by terrestrial microorganisms. 
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APPENDIX 


Consider a plane-parallel slab of total optical thickness 7, and let I(r, u) be the intensity of radiation at depth 7 
with p the cosine of the angle between the propagation direction and the forward normal. Assume isotropic 
scattering and let õ be the single-scattering albedo, which may be a function of 7. The equation of radiative 
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transfer (Chandrasekhar 1960; Irvine 1975) then reads 


d l L ~ a ¢yr Là ~ 
pw) = Tw) + Gok] de'T'(, W) + Bo exp (~m) (AD 
-1 


for an incident beam of unit intensity with direction cosine po, where 7’ is the intensity of the diffuse radiation and 
exp (— 7/0) that of the attenuated incident beam. The Eddington approximation can be obtained by substituting a 
function of a single variable 7,.'(r) for I’(7, u) when O < u < 1 and another function J_'(7) when —1 < p < 0. 
Let 


k = [3(1 — ão)” (A2) 
and let J’(7) be the average over all u of I(r, p). For constant x and a semi-infinite medium (rọ —> 00), one finds 


— 1+ QB, b = —> #0. 
FO 7? ET. 


The effective reflection coefficient R = 1 — f, where fis the fraction of the energy flux that is absorbed, is given by 


J'(r) = aexp (—«r) — bexp(—po™*7), (A3) 


PO i A4 

C+ Kpolll + Ci ao 

Generalizations for anisotropic scattering are discussed by Chandrasekhar (1960), Prinn (1970), and Irvine (1975). 

For isotropic scattering Prinn used equation (A3) for normal incidence (uo = 1) in his numerical work. For 
general uo but with x « 1, equation (A4) reduces to 


Suo) = 1 — Repo) X «(2/3 + uo). (A5) 


For many purposes one is interested only in results averaged over a whole daylight hemisphere, i.e., averaged 
over Ho from 0 to 1. Such averages are represented by the single problem we shall consider here, where the incident 
intensity [,,,(u) equals unity for all p between 0 and 1. In the spirit of the Eddington approximation we replace the 
problem by one for intensity (r, u) without explicitly isolating the attenuated incident beam: By analogy with 
equation (A1) we have 


pF Tew) =I) + ok | dI), (A6) 
which leads to the coupled equations 
FO JON- ad), ZO -= -Fo), (A7) 
where 
pl 1 1 
He) =4[ dle, FG) = duld, K= dI. (A8) 


We again make the approximation (but here for all the radiation) that I(r, p) = I-(r) for u > 0 and I(r, p) = 
I_(r) for u < 0, so that 


J=3K=%l, +1), Fai, -—I_). (A9) 


In this problem, any explicit consideration of the incident beam is simply replaced by the boundary condition that 
I,(7 = 0) = 1 and the effective reflection coefficient R = 1 — fis simply given by J_(r = 0). 

For general &,(7) and with x(7) defined by equation (A2), the use of the Eddington approximation (A9) in 
equation (A7) leads to the eikonal equation 


2 
E — eolo =0 (A10) 
and to the further relation 
W, +I)=3 s q- — L). (Ail) 
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We first solve these equations for the case of constant « for a slab of total thickness 7): We assume that the effective 
reflection coefficient Rp = I_(70)/I.(7o) is given as one boundary condition, the other boundary condition for the 
incident beam being J,(0) = 1. The general solution of the coupled equations (A10) and (A11) is 


I(r) = a{(3 + 2«) exp (—Kr) — (3 F 2x) exp [k(7 — 279) } , (A12) 
with a and y disposable constants. The requirement of 7_/I, = Ro at the second boundary (T = 79) leads to 


_ 8 — 2x) — Ro(3 + 2x) 
~ B= 2x) — R — 2k)’ 


and the condition 7,(0) = 1 then gives a. The final expression for the intensity J (averaged over direction cosines) 
is then 


(A13) 


_ {exp (—«r) — y exp [k(r — 270)]} 
he Rare ey ee 2k}y exp (— 2x79) at) 
Another quantity of interest is the effective reflection coefficient R, of the whole slab at the surface of incidence, 
which is oc J_(0), 


(3 — 2x) — y(3 + 2x) exp (— 2KTo) | 
i= (3 + 2x) — y6 — 2x) exp (—2K7>) (A19) 

Two limiting cases are of interest for the solutions given by equations (A12) to (A15). Consider first the case 
with Ro = 0 in the limit of x — 0 with 7o kept finite. In this case equations (A14) and (A15) reduce to 


4 _ 3r- 7) +2. 


a RT J(r) = 37, + 4 


(A16) 


The other case is the limiting case of To —> œ for any constant, nonzero x. Independently of the value of Ro, this 
case gives 


ha K < 1, the value of fin equation (A17) agrees with the average over 0 < po < 1 of equation (A5) which is 
= (4x/3). 

No exact solutions of the eikonal equation (A10) exist for general x(r), but if « varies sufficiently slowly, one is 
tempted to use the WKB approximation. For a semi-infinite medium, only one of the two WKB solutions survives, 
and we have 


F(a) œ [k(7)]-"2 exp |- Í ' dr'x(r’)| . (A18) 


Note that substituting a varying x(r) into equation (A17) (or into eq. [5] of Prinn 1970, which is the equivalent of 
our eq. [A3] for normal incidence instead of diffuse illumination) is not a good approximation to equation (A18) 
if x varies appreciably, even if the variation is slow. Furthermore, the requirement for the validity of the WKB 
approximation is that (dk/dr) « r?°, which is not satisfied by practical cases for «(r), such as the example in Figure 
2. An inelegant but reliable numerical method proceeds as follows: Subdivide the total slab into a large but finite 
number of thin slabs (with Ar % 0.3, say) and replace «(r) for each thin slab by its average value over this slab. 
Start with the Jast thin slab, assume Roe (3 — 2x)/(3 + 2x), and evaluate equations (A13) to (A15) with the slab 
thickness 7o replaced by Av. The value for R, obtained from equation (A15) for the last slab is then used as Ry for 
the penultimate thin slab, and J and R, are again evaluated, and so on iteratively until the first thin slab. The solid 
curve labeled J is the result of such a numerical solution for J(7), and the effective reflection coefficient R, for the 
whole medium (obtained at 7 = 0) is R, ~ 0.85. The dashed curve is the result for /(7) in a different problem—no 
absorption for 7 < T and complete absorption for r > 79, with rọ © 7.5 (the value of r for which «x = 0.5 in the 
actual problem). Equation (A16) gives the solution for J(7) for this problem and also yields R; ~ 0.85. Equation 
(A16) for R; should be a fairly good solution for any «(7), with ro the value where « ~ 0.5, as long as To > 1 and 
« is small for r appreciably less than 7». 

The accuracy of the numerical solution described above is limited (apart from computational errors) by errors 
inherent in the Eddington approximation, including the use of a constant coefficient in equation (A11). For a 
single slab this approximation is excellent when « « 1, although other two-stream approximations are better 
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(Sagan and Pollack 1967; Irvine 1975) when &o « 1. For the problems with 7 > 1, discussed here, « is small near 
the incident surface and increases only slowly with 7 (at first). The results for J(7) and R, are most sensitive to 
layers with moderately small 7 and are insensitive to inaccuracies in the Eddington approximation at large r. 
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